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ABSTRACT 

We present the first radial velocity curve of the companion star to BW Cir which demonstrates 
the presence of a black hole in this X-ray transient which recorded outbursts in 1987 and 1997 (and 
possibly 1971-2). We identify weak absorption features corresponding to a G0-5III donor star, strongly 
veiled by a residual accretion disc which contributes 61-65 % of the total light at A6300. The Doppler 
motions of these features trace an orbit of P = 2.54448 days (or its 1-yr alias of P = 2.56358 days) 
and velocity semi-amplitude K2 — 279 ± 5 km s""'^ (or K2 — 292 ± 5 km s~^). Both solutions are 
equally possible. The mass function implied by the shorter period solution is f {M) — 5.75 ± 0.30AfQ 
which, combined with the rotational broadening of the tidally locked companion {Vsini = 71 ± 4 
km s""'^), yields a compact object mass of Misin^i — 7.34 ± 0.46Mo. This is substantially above 
the mass of a neutron star under any standard equation of state of nuclear matter. The companion 
star is probably a G subgiant which has evolved off the main sequence in order to fill its Roche lobe. 
Remarkably, a distance of >27 kpc is inferred by the companion's luminosity and this is supported 
by the large observed systemic velocity (7 = 103 ± 4 km s~^) which requires such a distance in order 
to be consistent with the Galactic rotation curve. 

Subject headings: accretion, accretion disks - binaries: close - stars: individual: BW Cir (=GS 1354-64) 
- X-rays:stars 



1. INTRODUCTION 

Low mass X-ray binaries (LMXBs) are interacting bi- 
naries where a low mass star transfers matter onto a 
neutron star or black hole (e.g. Charles & Coe 2004). 
Mass transfer takes place through an accretion disc where 
angular momentum is removed and gravitational poten- 
tial energy converted into (high-energy) radiation. Ac- 
cretion discs evolve toward stationary states where the 
mass transfer rate through the disc, M, adjusts towards 
the value of M2, the mass transfer rate supplied by the 
donor star. M2 is driven by the binary/donor evolu- 
tion and if M2 < Merit ~ 10~^ Mq yr^^ then mass 
transfer instability cycles (outbursts) are triggered (van 
Paradijs 1996). There are -200 bright {hx ^ lO^^-lO^s 
erg s~^) LMXBs in the Galaxy and most of them har- 
bour neutron stars as implied by the detection of X- 
ray bursts/pulsations. On the other hand > 70 % of 
transient LMXBs are known to harbour accreting black 
holes, as demonstrated by dynamical studies of the faint 
companion stars, which is mainly possible when X-rays 
switch off (but see Hynes et al. 2003). We currently have 
dynamical evidence for 16 black hole LMXBs, with mass 
functions ranging from 0.22 ± 0.02 Mq up to 9.5 ± 3.0 
Mq (see Charles & Coe 2004, McClintock & Remillard 
2004) . The mass spectrum of black holes is of crucial as- 
trophysical significance in constraining supernovae mod- 
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els and the equation of state of nuclear matter. Clearly 
both more, and more accurate, black hole mass deter- 
minations are required before these issues can be fully 
addressed. 

BW Cir is the optical counterpart of the X-ray tran- 
sient GS 1354-64, discovered in 1987 by the Ginga satel- 
lite (Makino ct al. 1987). It displayed X-ray proper- 
ties reminiscent of BH transients i.e. a combination of a 
soft multi-BB component, with inner disc temperature of 
'-^^0.7 keV, plus a hard power-law tail with photon index 
2.1 (Kitamoto et al. 1990). BW Cir went into outburst 
again in 1997 but remained in a low/hard state through- 
out (Revnivtsev et al. 2000, Brocksopp et al. 2001). 
Interestingly, its sky position coincides with two older 
recorded X-ray transients, Cen X-2 and MX 1353-64, dis- 
covered in 1966 and 1971-2 respectively (see Kitamoto 
et al. 1990). The former was the first X-ray transient 
ever discovered, with a peak soft X-ray luminosity of ~ 8 
Crab. If this activity were attributed to the same source, 
then it would make BW Cir a BH transient with one of 
the shortest recurrence time (~ 8-10 years). However, 
the X-ray properties of these events were markedly dif- 
ferent and hence, if caused by the same source it would 
indicate that it has displayed at least four distinct X-ray 
states. 

BW Cir settles down in quiescence at R=20.5 (Martin 
1995) and as yet no indication of the orbital period exists. 
A possible period of ~ 46 hr and 0.3-0.4 V-mag ampli- 
tude was reported during the 1987 outburst (Ilovaisky et 
al. 1987) but it has not been confirmed. On the other 
hand, Martin (1995) suggests a very tentative 15.6 hr 
periodicity with ~0.1 R-mag amplitude from quiescent 
data although the folded light curve does not mimic the 
classic ellipsoidal modulation. Here we present the first 
spectroscopic detection of the companion star in BW Cir 
and the analysis of its radial velocity curve. This pro- 
vides the first dynamical probe for the presence of a black 
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Fig. 1. — pcriodogram of the radial velocities obtained by 
cross-correlation with the G5III template HD 62351. The best 
periods are at 0.393 day~^ (2.5445 d) and its 1-year alias at 0.390 
day-i (2.5641 d). 



hole in this historical transient X-ray binary. 

2. OBSERVATIONS AND DATA REDUCTION 

We have observed BW Cir using the F0RS2 Spectro- 
graph attached to the 8.2ni Yepun Telescope (UT4) at 
Observatorio Monte Paranal (ESO) on the nights of 22- 
23 June 2003, 14-15 and 25-27 May 2004. A total of 
55 spectra were collected with integration times varying 
between 1800 and 2000s (depending on atmospheric con- 
ditions). The R1200R holographic grating was employed 
which, combined with a 0.7 arcsec slit, produced a wave- 
length coverage of AA5745-7230 at 70 km s'^ (FWHM) 
resolution, as measured from gaussian fits to the arc lines. 
A He-|-Ne-|-Hg-|-Cd comparison lamp image was obtained 
with the telescope in park position to provide the wave- 
length calibration scale. This was derived by a 4th-order 
polynomial fit to 37 lines, resulting in a dispersion of 
0.74 A pix^^ and an rms scatter of 0.04 A. Instrumen- 
tal flexure in our target spectra was monitored through 
cross-correlation between sky spectra and it was always 
less than 33 km s~^. These velocity drifts were removed 
from each individual spectrum, and the zero point of the 
final wavelength scale was fixed to the strong 01 sky line 
at A6300.304. Two additional 1800s spectra were ob- 
tained on the nights of 23-24 Aug 2000. Here we used 
the R600R grism in combination with a 0.7 arcsec slit 
which resulted in a resolution a factor of two lower. 

We also observed several radial velocity and spectral 
type standards with exactly the same instrumental con- 
figuration on the nights of 17 May 2003 and 27-28 May 
2004. These cover spectral types F5III-K7III. All the im- 
ages were de-biased and flat-fielded, and one-dimensional 
spectra extracted using conventional optimal extraction 
techniques in order to optimize the signal-to-noise ratio 
of the output (Horne 1986). 

3. THE RADIAL VELOCITY CURVE AND SYSTEM 
PARAMETERS 

The 57 individual spectra were prepared for cross- 
correlation analysis by subtracting a low order spline fit 



to the continuum and rebinning the wavelength scale into 
a constant step of 34 km s^^ per pixel. Radial velocity 
points were extracted from our individual target spectra 
by cross-correlation with the radial velocity templates in 
the spectral range AA5950-7200, after masking out all the 
emission, telluric and IS absorption lines. The radial ve- 
locity points show clear night-to-night velocity changes 
which strongly points towards a 2d orbit. This was 
best seen on the nights of 14-15 May 2004, when two con- 
tinuous 9hr runs could be obtained. In order to search 
for the orbital period we have performed a power spec- 
trum analysis on the radial velocities obtained with the 
G5III template and the results are displayed in Figure 
1. Frequencies longer than 0.5 cycles d~^ (periods < 2 
days) can be ruled out since they provide xt > 15. The 
minimum xt (=2.98) is found at 0.3930 d~^, correspond- 
ing to P = 2.54446(6) days. However, we cannot rule out 
the 1-year ahas of P = 2.56358(7) days with xl = 3.01. 
All other aliases have > 5 and can be rejected since 
they are not significant at the 99.99 %. 

We have performed least-squared sine-wave fits to the 
radial velocity curves obtained with all templates and 
find that spectral types in the range F5-G8III give the 
best fits with xt in the range 2-3. We decided to adopt 
the fitting parameters of the G5III template since this 
provides the closest representation to the observed spec- 
trum (see below) i.e. 

7 = 103.1 ± 4.2 kms"^ P = 2.54448 ± 0.00015 d; Tq = 
2453140.985 ± 0.013; K2 = 279.3 ± 4.9 kms"^ or 

7 = 94.8 ±4.2 kms"^;P = 2.56358 ± 0.00015 d;To = 
2453140.939 ± 0.013; K2 = 291.7 ± 5.3 kms"^ 
where Tq corresponds to standard phase i.e. the in- 
ferior conjunction of the optical star. All quoted errors 
are 68 % confidence and we have rescaled the errors so 
that the minimum reduced x^ is 1-0. The 7- velocity has 
been corrected from the template's radial velocity and 
we note that it is unusually large for black-hole binaries, 
and only comparable to the long period microquasars 
GRO J1655-40 and V4641 Sgr. Figure 2 presents the ra- 
dial velocity points, folded on the 2.54 day period, with 
the best sine fit superimposed. The mass function of 
BW Cir. for the case of the short period solution, is 
/(M) = Misin-^/(l + g)2 = PKl/2i:G ^ 5.75±0.30Mq 
, which defines a lower limit to the compact object's 
mass. This is substantially above any neutron star mass 
defined by all standard equations of state of condensed 
matter and hence we conclude that BW Cir contains a 
black hole (Rhodes & Ruffini 1974). 

Furthermore, our high instrumental resolution (70 km 
s~^) enables us to measure the rotational broadening of 
the companion star and we can use this information to 
set a more stringent limit to the black hole mass. We 
have broadened our G5III template from 50 to 100 km 
s~^ in steps of 5 km s~^, using a Gray profile (Gray 
1992) with a linear limb-darkening law with coefficient 
e=0.62, interpolated for A=6300 A and Toff ~ 5000 K 
(see Al-Naimy 1978). The broadened versions of the 
template star were multiplied by fractions / < 1, to ac- 
count for the fractional contribution to the total light, 
and subsequently subtracted from the grand sum spec- 
trum of BW Cir. The latter was produced after aver- 
aging the 55 high resolution spectra of BW Cir in the 
rest frame of the companion using the orbital solution 
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Fig. 2.— Radial velocity points of BW Cir folded on P = 2.54448 
days and Tq = 2453140.985 with the best fitting sinusoid. 



above. A test on the residuals yields V smi = 71 ± 4 
km s~^ (we have repeated the same experiment with all 
the templates and Vsini always ranges between 59-71 
km s~^, see table 1). In the case of a tidally locked 
Roche-lobe filling star, Vsini relates to the _R'2-vclocity 
and the mass ratio q — M2/M1 through the expression 
Vsini ~ 0.462i^2g^/^(l + (Wade & Home 1988), 

from which we obtain q = 0.13 ± 0.02 and, therefore. 
Ml sin^ i = 7.34 ± O.46M0. Since BW Cir does not ex- 
hibit X-ray eclipses, the inclination angle must be i < 77° 
(for q ~ 0.13) and, consequently, the masses of the black 
hole and the companion star are Mi > 7.83 ± O.SOMq 
and M2 > 1.02 ± 0.17Mq respectively. 

For the longer period solution, the black hole case is 
even stronger, with /(M) = 6.60 ± O.36M0, q = 0.12 ± 
0.02, Ml sin^i = 8.28 ± O.54M0, Mi > 8.95 ± O.58M0 
andM2 > l.O7±O.19M0. Both the long period and large 
mass function are reminiscent of the black hole transient 
V404 Cyg (Casares, Charles & Naylor 1992). 



4. THE NATURE OF THE COMPANION STAR AND 
DISTANCE TO BW CIR 

In order to refine the spectral classification we have 
averaged our individual spectra in the rest frame of the 
companion star, using our first orbital solution, and com- 
pared it with our optimally broadened spectral type tem- 
plates in the range F5-K7III using a minimization 
routine (see Marsh, Robinson & Wood 1996 for details). 
The minimum is obtained for spectral types GOIII- 
G5III, which contribute 39-35 % to the optical flux in 
our spectral range (see table 1). Figure 3 presents the 
averaged spectrum of BW Cir in the companion's rest 
frame, together with our G5III template^, broadened by 
Vsini = 71 km s~^. 

A normal G0III-G5III typically has M2 = 2.1 - 2.4M0 
and i? ~ 5 — 8i?0 (Gray 1992). However, the Roche lobe 
equivalent radius of a 2.1—2AMq star in a 2.5d orbit is ~ 

^ HD 62351 is a close binary with a separation of 0.4" , resolvable 
through speckle. It is classified in the V-band as a G5III or G5/6IV 
(Christy & Walker 1969, Houk & Smith-Moore 1988). We confirm 
this classification in the R-band based on the line ratio Ha/ X6498, 
which is very sensitive to T^g (see table 1). 
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Fig. 3. — Doppler corrected average of the 57 individual BW 
Cir spectra in the rest frame of the companion star and the G5III 
template, broadened by 71 km s^^. Note the extreme weakness of 
the metallic features from the companion star, as compared to the 
template. 



4.6 — 4.8J?0 and hence the companion is probably a lower 
mass, but more evolved star. Since T^jf is constrained 
by the spectral classification (~ 5000 K), and the radius 
of the companion by the Roche lobe geometry'' 
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we can apply the Stefan-Boltzmann relation to constrain 
the distance to BW Cir. Our orbital solutions yield 
M2 > 1.02 M0 and hence i?2 > 3.6 R0. Therefore the 
Stefan-Boltzmann relation gives M^oi < 2.6 which, com- 
bined with the bolometric correction and (V-R) colour 
for a G5III star (Gray 1992), gives Mr < 2.11. The 
dereddened magnitude of the companion star can be esti- 
mated from the observed quiescent magnitude (R=20.5), 
corrected for reddening {Eb-v ^ 1; Kitamoto et al. 
1990) and our determined 65 % veiling. These yield 
R ~ 19.3. Therefore, the distance modulus relation pro- 
vides a lower limit to the distance of 27 kpc which is 
substantially larger than the 10 kpc value proposed by 
Kitamoto et al. (1990). Even so, we note that this makes 
Lx(l - lOfceV) = 3.5 X l^f^ {d/lQkpcf > 2.6 x lO^^ erg 
s~^ at the peak of the 1987 outburst, still confidently 
below the Eddington luminosity for a 10 M© black hole. 

The distance has important implications for the ob- 
served systemic velocity. Wc have calculated the ex- 
pected radial velocities of BW Cir, due to Galactic dif- 
ferential rotation, for a distance of 10 kpc and 27 kpc 
and these are -6 and 104 km s""'^ respectively (Nakanishi 
& Sofue 2003). The large distance scenario enables us 
to explain the observed 7-velocity without invoking any 
kick during the formation of the black hole. Although 
this scenario is the most attractive, we note the uncer- 
tainties involved in the IS extinction value and this re- 
sult requires careful consideration. Our spectrum shows 

^ This results from combining Kepler's Third Law with 
Paczynski's relation (Paczynski 1971). 
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Fig. 4. — H-R diagram with evolutionary tracks of solar metal- 
licity stars of masses 1.0-2.0 Mq, from Cassini, Salaris &; Castelli 
(2004). Dotted lines mark the space region allowed by our spec- 
tral type classification (T^ff = 5091 — 5743 K) and the Stefan- 
Boltzmann relation for the case of a Roche lobe filling star with 
M2 > 1.02 M0 and P=2.544 days. 



several IS absorption features which can be used to ob- 
tain an independent estimate of the reddening. The 
strongest ones are the Nal doublet for which we mea- 
sure EW(D2)=1.37± 0.04 A and EW(D1)=1.66±0.04 A. 
The empirical calibration of Bardon et al. (1990) yields 
Eb-v = 0.25 X EW = 0.78 but this should be regarded 
as a lower limit because NaI(D2) is somewhat diluted by 
the broad Hel A5875 emission from the accretion disc. 
We have also measured the EW of the diffuse IS band at 
A6203 and find 0.24 ±0.02 A. Using the empirical cali- 
bration between EW and reddening from Herbig (1975) 
we get Eb-v — 0.84 ± 0.10. These numbers are con- 
sistent with Kitamoto et al.'s reddening of Eb-v ~ 1, 
which was estimated by comparing the observed outburst 
colour of the optical counterpart with typical values of 
luminous LMXB given in van Paradijs (1983). This as- 
sumed that the X-ray irradiated disc dominated the op- 
tical emission, which is confirmed by our results here and 
hence give further support to our estimated distance of 
>27 kpc. We note that this distance argues against BW 

* Pietriferni e t al. 2004, in preparation; 

p:tp: / /www.te. astro. it /B ASTI| 



Cir and Cen X-2 being the same object since it would 
make the 1966 outburst substantially super-Eddington 
(by more than an order of magnitude) for a 10-20 Mq 
black hole. However, it does make BW Cir similar to the 
black hole LMXB GX339-4, where analysis of the Nal IS 
aborption features at high resolution has revealed com- 
plex velocity structures, consistent with a large distance 
of > 15 kpc (Hynes et al. 2004). We also note that this 
implies Lx ^ 0.1 LEdd for the 1987 hard X-ray outburst 
which makes the most luminous hard state yet seen in a 
BH transient. 

Additional constraints are provided by the position of 
the companion star in the H-R diagram, as compared 
with evolutionary sequences of single stars. Our spec- 
tral analysis indicates that the companion is well fitted 
by a G0-5III template and hence T^// ~ 5091 - 5743 
K (Gray 1992). On the other hand, combining eq. 1 

with the Stefan-Boltzmann relation gives log ^ 
4 X logTe// - 13.92, where we have used M2 > 1.02 Mq. 
This defines a region in the H-R diagram (see Figure 4) to 
compare with evolutionary tracks of 1-2 Mq stars, com- 
puted using the code BASTI* for the case of Population 
I stars of solar metallicity. The diagram shows that the 
companion star is located in the Hertzsprung gap or the 
lower part of the giant branch for M2 > 1.6—1.8 Mq. 
Since the luminosity of a Roche lobe filling star depends 
weakly on M2 , the companion cannot be much more mas- 
sive than ~ 1.6Mq (for instance logL2 = 1-06 Lg for 
M2 = 2.0AfQ). If it were a 1.6 Mq star with Teff - 5000 
K, then logL2 — 1.0 Lq and the distance would be ~ 33 
kpc. However, we note the limitations of this crude es- 
timate because the donor stars in X-ray binaries are not 
in thermal equilibrium and their H-R tracks may deviate 
significantly from single star evolution. 
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TABLE 1 

Radial Velocity solutions and Spectral Type classification 



Template 


Spect. Type 


//a/A6498 


7 


P 


To 


K2 


V smi 


f 










(km s-l) 


(days) 


(+2453140.0) 


(km s-l) 


(km s-i) 






HD 182901 


Fsni 


8.2(1) 


103(3) 


2.54447(12) 


0.987(9) 


276.4(3.4) 


66(4) 


0.47(2) 


2.50 


HD 83273 


GOHI 


6.0(1) 


102(3) 


2.54445(9) 


0.983(8) 


279.3(3.2) 


69(4) 


0.39(1) 


2.47 


HD 62351 


G5HI 


3.7(1) 


103(3) 


2.54448(9) 


0.985(8) 


279.2(2.9) 


71(4) 


0.35(1) 


2.42 


HR 6915 


G8HI 


2.7(4) 


104(2) 


2.54451(8) 


0.990(7) 


278.0(2.7) 


71(4) 


0.25(1) 


2.56 


HR 6925 


K3HI 


1.6(1) 


109(2) 


2.54451(7) 


1.000(6) 


279.2(2.0) 


59(4) 


0.15(1) 


2.58 


HD 184722 


K7HI 


1.3(1) 


110(2) 


2.54448(7) 


0.990(6) 


283.7(2.1) 


61(4) 


0.13(1) 


2.66 
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